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Abstract 

We generalize linearized Gouy-Chapman-Stern theory of electric double layer 
for nanostructured and morphologically disordered electrodes. Equation for 
capacitance is obtained using linear Gouy-Chapman (GC) or Debye-Hiickel 
equation for potential near complex electrode/electrolyte interface. The ef- 
fect of surface morphology of an electrode on electric double layer (EDL) 
is obtained using "multiple scattering formalism" in surface curvature. The 
result for capacitance is expressed in terms of the ratio of Gouy screening 
length and the local principal radii of curvature of surface. We also include 
a contribution of compact layer, which is significant in overall prediction of 
capacitance. Our general results are analyzed in details for two special mor- 
phologies of electrodes, i.e. "nanoporous membrane" and "forest of nanopil- 
lars". Variations of local shapes and global size variations due to residual 
randomness in morphology are accounted as curvature fluctuations over a 
reference shape element. Particularly, the theory shows that the presence of 
geometrical fluctuations in porous systems causes enhanced dependence of 
capacitance on mean pore sizes and suppresses the magnitude of capacitance. 
Theory emphasizes a strong influence of overall morphology and its disorder 
on capacitance. Finally, our predictions are in reasonable agreement with 
recent experimental measurements on supercapacitive mesoporous systems. 
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1. Introduction 

Electrochemical study on curved nanostructured surfaces (CNS) is due 
to the requirement to develop efficient energy generating and storage de- 
vices [lj, electro- mechanical systems j^] and their applications in nanofluidics 
[3| . Recently, the electrochemical supercapacitive behavior of porous carbon 
materials [If, e.g. carbide derived carbons (CDC) (jlS, activated carbon 
(ACs), graphitic carbon, carbon nanotubes (CNTs)j2| etc., have acquired 
intense focus. These supercapacitive systems mainly based on the electrical 
double layer (EDL), requires surfaces with high specific area and volume with 
proper pore size and shape control for the efficient access of ions to obtain 
high energy storage as well as high power. Hence, the electrical double layer 
(EDL) formed at the complex nanostructured or disordered interface is a 
major focus for the research^, |, @, 0, @, S, M, E2, Q El . 



Recently, Gogotsi and coworker were successful in developing supercapac 
itor with well controlled pore sizes in porous carbide derived carbon material 
(CDC) Q, i [l~4| which shows an anomalous behavior in the capacitance (lil 



These experimental results show three regimes in the capacitance vs pore size 
data[H, a, H] : (i) shows a nonlinear increase in capacitance, (ii) shows the 
transition from micropores to mesopores capacitance with a minimum, (iii) 
shows the anomalous increase and a maximum in capacitance of pores below 



1 nm. Huang and coworker [16|, LL7J, ll8| proposed a heuristic model for the 



capacitance of such a problem with an assumption of cylindrical pores. They 
included the effect of finite pore sizes and proposed three models for three 
different pore size regimes: (i) micropore regime (< 2nm)- the electric wire- 
in-cylinder capacitors (EWCCs) model, [13] (ii) mesopores regime ( 2-50 nm)- 
the electric double cylinder capacitors (EDCCs) model [l7| and for pores in 
macropore regime (> 50 nm) where curvature is not longer significant, the 
parallel plate capacitor model was proposed. The capacitance (C) per unit 



area (A) of EDCCs model is given as 17 



and for EWCCs is 11 



C,A ~ b Hb/(b- d)\ (1) 



C,A ~ b \n(b/a ) (2) 



2 



in which A is the area of cylindrical pore, e r is the electrolyte dielectric con- 
stant, eo is the permittivity of vacuum, d is the effective thickness of the 
electric double layer (the Debye or Gouy length), do is an effective size of the 
counter-ions in core and b radius of the outer cylinder. Computer simulations 
are becoming standard tools to study EDL on planar geometri es | l9l . |20| . u 
but recently applied to curved and porous systems 22j, [23j, [24|, [25] , with ac- 



counting ion size, ion- ion correlation [20L |2l| and ion-solvent interaction [26 



27j. Also simulations of EDL with accounting the dependency of the elec- 



trolyte dielectric permittivity on the local electric fields [28( in sphere and with 



electrodes made of closely packed monodispersed mesoporous spheres [29 



These simulations do not account for the general geometrical features and 
morphological fluctuation in the electrode. A general theory must account 
for the geometrical and morphological feature of an electrode that affects the 
capacitance. 

The porous electrode materials are ubiquitously complex. The modelling 
such systems is difficult, and complexity arise due to pore structure with inter 
connected three-dimensional connectivity of pores of nonuniform shape and 
size. The morphology of the pores and particulate materials may be idealized 
to various forms - prolate, oblate, ellipsoidal, spheroidal, tubular, etc. These 
forms and morphologies can easily be identified with their local curvatures. 
But for the real pore space will be a combination of these idealised shapes 
or can be looked upon as fluctuations in pore structure around one of these 
idealized shapes hence identified with statistical property of their curvatures. 
The electrochemical capacitance is strongly influenced by the morphology of 
nanoporous electrode materials 3(|. Hence, there is a need for systematic 
rigorous theory to understand the capacitance of nanostructured and disor- 
dered electrodes [3l|] which interpret capacitance behavior of such electrodes 



within the frame work of electric double layer theories |32l. 133 



In this article, we develop (linearized) Gouy- Chapmann- Stern level theory 
for an arbitrary morphology. Simplicity of the model is maintained through 
linearization while ignoring the microscopic interaction, viz, ion-solvent in- 
teraction and ion-ion correlation, which may influence the ion- distribution. 
The theory is based on segregation of the compact layer and diffuse layer 
regions with an assumption of validity of the linearized Poisson-Boltzmann 
equation in the diffuse layer region while compact layer corrections are in- 
cluded at the level of local capacitance density. In section 2, we understand 
the effect of geometry and topology on the diffuse layer capacitance. In 
section 3 we have included the compact layer contribution to diffuse layer ca- 
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pacitance on an arbitrary curved surface electrode. In section 4 we develop 
a theory for the nanostructured electrode with random morphology, viz, cur- 
vature fluctuation in a porous system arises due to intrapore roughness and 
interpore size variations. General results are applied for detailed analysis for 
the porous membrane and forest of nanorods electrode. Comparison with 
recent experimental data of capacitance to pore size is presented in section 
5. Finally, the summary and conclusions are discussed in section 6. 




Figure 1: Schematic model of EDL formed in arbitrary nanoporous electrode showing the 
Helmholtz layer and the diffuse layer. 



2. Model of Diffuse Layer near Curved Surfaces 

Historically, the concept of EDL is dated back to 1879. It was Helmholtz 
34J |. who first conceived the idea that adsorption of opposite charge ions 
from solution on an electrode surface with excess or deficiency of charge 
results in a situation like a parallel plate capacitor of opposite charge with 



few nanometer thickness. The capacitance density of Helmholtz layer |8[ is 

c H = (3) 

r H 

where eo is the dielectric constant of free space, e# is the HL relative dielectric 
constant, and r# the thickness up to the outer Helmholtz plane (OHP). Nagy 
et al. showed that the dielectric constant in HL capacitance (e#) should 
be the average of inner layer and diffuse layer dielectric constants due to 
the polarization charges induced on the boundary of the Helmholtz and the 
diffuse layer 3 • 
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Later Gouy[35l| and Chapman [36[ included the fact that the ions are mo- 
bile in the electrolyte solution due to thermal motions and developed a math- 
ematical treatment based on combined application of Boltzmann's energy 
distribution equation and Poisson's equation. However, the Gouy-Chapman 
model is known to over estimates of the EDL capacitance for planar elec- 



trode. The divergence was overcome by Stern [37J in 1924 where he combined 
in series both the HL capacitance (c# ) and the Gouy-Chapman capacitance 
of a planar electrode (cqc) as: ^~ — ^ + where cqc = ^g^cs/Ait and 
K c ff = k cosh(e/30/2). The Debye-Hiickel screening length, kT 1 — Id — 
(eksT /Ane^n^zf) 1 / 2 where T is temperature, ks the Boltzmann constant, 
e the dielectric constant of bulk electrolyte, is the number density of i-th 
ion with charge z«e and e is the electronic charge. 

Here we will like to develop a model for diffuse layer near a curved surface. 
Using Gauss law, the local charge density (a) in electric double layer is: 
a = (e/47t)(d<j)/dn); where <fi is the potential relative to the bulk solution 
and d/dn = n ■ V is the normal derivative (with unit normal vector (n) to 
the surface). The specific differential capacitance of a diffuse double layer 
is defined as cg = da/d(p , where O is the potential difference applied at 
the interface. Differentiating a gives capacitance density in terms of the 
electrostatic potential as 

c G = ±A. ( *) (4) 



An d(j)Q \ dn r 

here e is the dielectric constant of the bulk electrolyte and d<ft/dn = h ■ V</> 
the inward normal derivative of potential to the surface which is a functional 
of surface potential. 

The potential in Eq. H] can be obtained using Poisson-Boltzmann equa- 



tion (PBE) and Tessier and Slater 38j have shown the validity of PBE when 



the system size is reduced to nanoscopic dimensions. Daikhin et al. [39 
showed that result for the nonlinear PB theory of capacitance at a rough 
electrode electrolyte interface is related to result for the linear version of PB 
Equation. On similar lines to include electrodes of arbitrary morphologies, 
we developed a theory based on simplified linearized Poisson-Boltzman equa- 
tion. The effective nonlinear contribution of PB equation may be included in 
the linear solution by replacing the inverse Debye length (ac) with a potential 
dependent inverse effective Debye or Gouy length K eS = k cosh(e0/2fc B T) 



4G . |41| . The solution of Poisson's equation for the potential 4> relative to 
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the bulk electrolytes can be obtained under the condition of strong or weak 
potential approximations. Hence, the simplifying assumption used in our 
work is that the regions near an electrode are divided into a strong electric 
field region, i.e. the compact layer, and the region beyond the compact layer 
has a relatively weak electric field, i.e. the diffuse layer (hence assumed to 
obey linearized Gouy-Chapman equation or Debye-Hiickel equation ). 

The linearized PBE or the linearized Gouy-Chapman equation, for po- 
tential (0(r)) relative to the bulk solution, is written as 



(V 2 - K 2 )0(r) = 0. 



(5) 



where the interfacial potential w.r.t. bulk is taken as constant and the bulk 
electrolyte potential is assumed to be 0. The solution of the Helmholtz equa- 
tion or Debye-Hiickel equation for arbitrary geometries and boundary con- 
ditions (Dirichlet, Neumann or Robin) are known in the form of "multiple 



scattering expansions" [42|, |43|, |44|, |45|, |46[ . Also " multiple scattering expan 



sion in curvature" results accounting for the arbitrary geometry of the surface 

Adapting a sim- 



is analysed for the diffusion problems |46|, |47|, |48|, |49|, [5C| [51| 
ilar methodology to solve the electrostatic potential of electric double layer, 
we will obtain the capacitance of EDL near surfaces of arbitrary shapes in 
Debye-Hiickel or linearized Gouy-Chapman regime. The convergence of the 
screening length power series is applicable in a strong screening regime, viz, 
the Debye-Hiickel screening length, k -1 is smaller than any scale of radius of 
curvature. Through local and global curvatures brings in the geometric and 
topological features of surface are expected to play a role in the formation 
of the electric double layer near a curved surface. The expansions will be 
obtained in powers of k _1 / R, where R is a typical radius of curvature. 

We use the method of Green function in order to obtain the various order 
of scattering terms depending on surface curvature, the detailed calculation 
is shown in appendix. Now the capacitance density at ct-th point of surface 
in terms of Green function is obtained using Eq. H] as 



c G (a) 



' d 3 r> 



€K 



47T 



3 „,0 G(q + ,r') 
dn a 

dG (a,r') 



v 



dn r 



22 J dG (a,P) dG {/3,r') dS 



dn r 
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3 / dGo{cL,P)dG Q {M) dGpfrS, 

+i I ~ k k aapdbj — ■ 

on a orifs an 1 



(6) 



The local shape of the interface is given by mean (H a ) and Gaussian curva- 
tures (K a ) at the a-th point on the surface. Expressing the various orders of 
scattering terms are obtained in terms of two curvatures of the surface, the 
capacitance density of diffuse layer is (see appendix for details) 



c G (H a , K a ) = — 



l--H a -±- 2 {Hl-K a ) + 
K 2k z 



(7) 



where the local quantities H a and K a are defined as (dropping subscript for 
simplicity), H = (l/Ri + l/-R 2 )/2 and K = l/RiR 2 . Retaining Eq. [TJupto 
second order in curvature expansion where two curvatures satisfy: (H 2 — 
K) = 1/4 {l/Rx - l/R 2 f and (H 2 - K) > 0. Eq. clearly shows that the 
geometric dependence of the capacitance is controlled by second and third 
term through H and K. If H = and K = 0, the diffuse layer capacitance 
simplifies to the flat surface ( for small applied potential) Gouy-Chapman ca- 



pacitance |41] . The local capacitance result for an arbitrary curved electrode 
has three terms. The first term is dependent on the solution properties, sec- 
ond term is dependent on geometry (independent of ionic concentration) but 
third term represents the coupling between geometrical and solution prop- 
erties. Eq. [7] for diffuse layer capacitance can be interpreted as an effective 
parallel curved surface capacitor, where two curved surfaces are separated by 
an effective distance n~ l / [l — \H a — -^(H 2 — K a )\ . The local separation 
between electrode and effective surface in solution depends on the composi- 
tion as well the local curvature of the electrode. 

We illustrate use of Eq. [7] through its application to a cylindrical pore and 
a rod geometries. Fig. 2 shows the plot for the capacitance of these simple 
geometries and their dependence on size (2 r). Two curvatures for a cylindri- 
cal pore, mean curvature H = l/2r is constant and the Gaussian curvature 
K = 0. Similarly, for a cylindrical rod, H = — l/2r while K = 0. In Fig 
2 we show two cases of capacitance, viz. linearized Gouy-Chapman (LGC) 
model capacitance for curved surfaces (dash lines) and LGC capacitance after 
Helmholtz layer size correction we call this attuned LGC capacitance (AGC) 
for curved surfaces (solid lines), viz, for a cylindrical pore the attuned curva- 
ture is H = 1/(2 (r — th)) and for a cylindrical rod the attuned curvature is 
H = —\/{2{r+r H )). The LGC and AGC capacitance curves for the pore have 
smaller values compared to flat surface. In the case of rod (red lines) we see 
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2r(nm) 

Figure 2: Theoretical capacitance density vs diameter plots for cylindrical pore and rod 
electrodes. The linearized Gouy-Chapman (LGC) capacitance curved surface (dotted 
lines) and attuned LGC capacitance (solid lines) after inclusion of Helmholtz layer size 
correction for pore (red line) and rod (blue line) geometries. The plots are generated using 
rn = 0.35 nm, e H = Q and e = 38 for an organic electrolyte at 298 K. 



that the LGC capacitance value increases with reduced rod diameter while 
AGC this rise weakens due to larger attuned rod diameter. This finding for 
LGC is similar to Compton and coworkers report of enhancement of capaci- 



tance in diffuse layer for hemispherical and cylindrical electrode |32|, |52j. But 



with increasing sizes, all electrodes reach curvature independent capacitance 
value and merge with the planar capacitance limit. The point particle nature 
of electrolytes beyond the compact layer, allows it to be extended up to pore 
size equal to mono- layer of ions or say upto the compact layer size. In short, 
it is important to emphasize that our theory will be applicable to a porous 
system for pore sizes 2r > 2r# + Id, i.e. sufficient to accommodate compact 
as well as the diffuse layer inside them. Extension of curves below the pore 
diameter 2th + Id is allowed due to point particle nature in diffused layer, but 
in reality we will require contributions from volume exclusion, material space 
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charge corrections etc. for an accurate prediction of capacitance. Here we 
show the emergence of morphological influence on capacitance is mainly due 
to diffuse layer. It will be significantly seen in cases when the magnitude of 
the smallest radius of curvature (R) follows: m < R< {^Th + Id)- The exis- 
tence of minimum in the capacitance of highly confining pores is reported in 
literatures [H, [l4, 15, 17] but its explanation will require inclusion of material 



capacitance of the electrode. The large radius of curvature limit of capaci- 
tance of pore as well as rod has curvature independent capacitance regime 



(as of planar surface). The experimental observations [15|, l53( of maximum 
capacitance at a critical pore diameter which corresponds to the effective 
thickness of HL may require inclusion of space charge contribution form the 
electrode material in our theory. 

The diffuse layer capacitance given by Eq. may be represented for 
simple curved surfaces as 

c G = co(l + a^ + & !) (8) 

where Co = ^dI^d is planar capacitance density, (a, b) are pairs of dimen- 
sionless numbers depending on the shape of a nanoelectrode, viz. (0,0) for a 
plane, (-1/2, -1/8) for a tube, (1/2, -1/8) for a rod, (-1,0) for a cavity and 
(1,0) for a sphere. Effective attuned radius due to compact layer (r a = r±r#) 
on curved surfaces, where + and — signs are assigned to convex (e.g. sphere 
and rod) and concave (e.g. tube and cavity) surfaces, respectively. 

The diffuse layer capacitance obtained in Eq. [5] can be looked upon as 
the capacitance between two curved surfaces, e.g. electrode surface and an- 
other effective surface in solution with same symmetry at distance /d/(1 + 
a — + b -Tj-) (in the normal direction). More specifically, the capacitance of 
a cylindrical pore in Eq. [8] can be approximated as an effective concentric 
cylinder capacitor with their effective separation l D /(l — ^ — ^) and for a 

cylindrical rod electrode in presence of diffuse layer is equivalent to wire in 

; l 2 

a cylinder capacitor with their separation /_d/(1 + — g^r). 

Another important quantity of physical interest is the total diffuse layer 
capacitance of a nanostructured electrode, Co- This is obtained by integrat- 
ing the local diffuse layer capacitance density cq in Eq. [TJ with respect to 
the interfacial surface as 



C L> = l dSc G = ^- 



(9) 
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K gives the intrinsic property and topological information for a surface 
through theorem of Gauss-Bonnet j54| as ~K = J s dS/(RiR 2 ) = Am (N c - N h ) 
where N c is number of connected surface and Nh is the number of handles 
(genus or holes) in the surface. H gives the extrinsic property and knowledge 
of space. The morphological quantities are the geometric area A = J s dS, 
integral of mean curvature H = J s HdS, integral of square of mean curva- 
ture H 2 = J s H 2 dS and topological quantity, integral of Gaussian curvature 
K = J s KdS can be calculated both for regular and random geometries. 

3. Inclusion of Compact Layer Contributions in Local Diffuse Layer 
Capacitance 

The capacitance of the EDL at Gouy, Chapman and Stern level is usually 
obtained using an analogy of capacitors in series. Here we assume validity of 
ansatz of combining two components of capacitances in series. Hence local 
capacitance density (c s ) for curved interface is the local compact (c#) and 
local diffuse layer (cq) capacitance density in series and written as: 

± = i- + -L (10) 

C s C H C G 

Two limiting situations for the local capacitance are observed due to different 
functional dependence of the charge on the electrolyte concentration and 
electrode potential. 

Two components of capacities are estimated in our analysis by dividing 
the regions near an electrode into a strong electric field region- the capaci- 
tance estimated using the Helmholtz formula. The region beyond the com- 
pact layer is the diffuse layer and has a relatively weak electric field. The dif- 
fuse layer capacitance is estimated through AGC model for the curved inter- 
face. This approximation provides a workable simplicity to our model. Com- 
bining Eqs. [7] and [10] we obtained the attuned (linearized) Gouy-Chapman- 
Stern (AGCS) capacitance density for an arbitrary geometry electrode as a 
function of mean and Gaussian curvatures: 

niJ r K , c H c [l-Hl D -(H 2 -K)(ll/2)} 

CMK) = c H + c [l-Hl D -(H 2 -K)(ll/2)} (U) 

where cq = ck/Atc and ch = €H£o/ r H, where e# is the electrolyte dielectric 
constant in Helmholtz layer (HL) and eo is the permittivity of vacuum, rn 
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is the effective thickness of Helmholtz or compact layer. Here H and K are 
obtained at the Helmholtz surface (by local shifting of electrode surface in 
normal direction by distance r#), respectively. The diffuse layer starts be- 
yond the compact layer hence the curvature attuned at the Helmholtz layer 
is used in our calculation. The effective or attuned radii of curvatures at the 
diffused layer are accounted by subtracting or adding r# from radii of cur- 
vatures depending on local convexity of surface. Combining compact layer 
capacitance with AGC capacitance partially removes point particle limita- 
tions as well as the divergence under R/Ib —> limit. The local AGCS 
capacitance in Eq. [TTJhas three terms in the square brackets. The first term 
in square bracket is dependent on the solution properties, second term is 
dependent on geometry through the local mean curvature but third term 
represents the coupling between geometrical and solution properties. Third 
term has mixed mean and Gaussian curvatures, viz H 2 — K, which is equal 
to for sphere hence measure of deviation from local asphericity and follows 
constraint: H 2 — K > 0. Hence local coupling between geometry and ionic 
solution emerges due to deviation in local asphericity of surface. 



4. Capacitance for Random Morphology Electrodes 

Randomness arises due to disorder in porous material can be of three 
types: (1) pore roughness can be looked upon as intra-pore shape fluctuations 
along the contour of the pore, (2) inter-pore size fluctuations over different 
pores and (3) fluctuations in length of pores. First two contributions of 
randomness can be accounted through fluctuations in mean curvature. Some 
examples are: carbide derived carbons (CDCs) have intrapore fluctuations, 
template carbon materials (TCM) may have both intrapore fluctuations and 
interpore pore-size fluctuations, activated carbons (ACs) and hierarchical 
porous graphitic carbon material (HPGC)[55| may have intra-pore curvature 
fluctuations, pore size fluctuations and pore length fluctuations. Fluctuations 



in length of pore affects, mostly the dynamic response of pores |55|, |40j. This 
case is not required in calculations as one is not looking in the dynamic aspect 
of this problem. 

The most appropriate method for characterizing the complex disordered 
structure are based on the statistical approach. Depending upon the nature 
of disorder in nanostructured or porous electrode one can have various dis- 
tributions of curvature. We assume in our model of an electrode with small 
random surface roughness around its reference geometry and is characterized 
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in terms of an ensemble averaged mean and Gaussian curvatures, viz, (H) 
and (K) . The surface roughness is characterized in terms of their deviations 
from reference curvature: H — (H) and K — (K) with an ensemble aver- 
age mean values: (H — (H)) = and (K — (K)) = but finite variance 
((H — (H)) 2 ) (fourth order term in principle curvature, ((K — (K)) 2 ), is as- 
sumed to be negligible). In order to calculate the capacitance of disordered 
systems, we take the statistical ensemble of various configurations in cur- 
vature space (for which the curvatures are distributed). Now the ensemble 
averaged capacitance density of diffuse layer is obtained from Eq. [7] as 



where (•) represent an ensemble average over all possible random surface con- 
figurations and Id = ft" 1 - Eq- H2]has three terms: first term is dependent 
on only solution ionic concentration, second term is purely dependant on en- 
semble averaged mean curvature while third term is due to coupling between 
ionic solution and electrode morphology. The mean capacitance density (ca) 
of diffuse EDL will be strongly affected whenever the radii of curvature will 
be comparable to diffuse layer thickness, which is a strong function of elec- 
trolytes concentration. 

Deviation in an ensemble average capacitance density is proportional to 
the difference of the curvature of surface and ensemble average curvature. 
The capacitance deviation of diffuse layer (5cg = cq — (cg)) is 



Sc G = c [l D (H - (H)) - (l D /2) (H 2 - (H 2 )) + (l 2 D /2) (K - (K))] (13) 



The average deviation of capacitance is zero, i.e. (Scq) = 0. The mean 
square deviation of difference of capacitance density is obtained from Eq. [13] 



An important feature of practical surfaces is the presence of small surface 
roughness characterized as small fluctuations in surface curvatures about its 
reference geometry. We have truncated higher-order curvature contributions 
in Eq. [TH It is obvious from Eq. [TH that the ensemble average of square 
of capacitance deviation is directly proportional to the ensemble average 
of square of deviation in mean curvature. Rewriting Eq. [14] in terms of 
a coefficient of morphological fluctuations or relative variance of the mean 



(cg) = co \l-l D (H)- l f((H 2 )-(K)) + ... 



(12) 



as 




(14) 
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curvature, 7 2 = ((H — (H)) 2 } / (H) 2 , we have 



<*4> = (£) 2 7W + o(tf 4 ) 



(15) 



where 7 is a coefficient of morphological fluctuations and is a measure of 
relative deviation from the mean curvature. This variance in capacitance 
density is independent of concentration of electrolyte. Coefficient 7 < 1 for 
weakly fluctuating interface and 7> 1 for strongly fluctuating interface. 

As mentioned before porous materials have complex spatial structures 
and are characterized by using morphological measures, viz, the mean (mi- 
croscopic) geometric area of an electrode is A, ensemble average mean curva- 
ture (H), ensemble average square of mean curvature (H 2 ), and the ensemble 
average of Gaussian curvature (K). For random ergodic fields [56| which are 
statistically homogeneous over various configurations, the (large) surface av- 
erage over morphological quantities, viz, H/A, H 2 /A and K/A (as defined 
earlier), can be related to the ensemble average morphological quantities as 



(H) = XmiA^H/A, (H 2 ) = linu_>oo# 2 M and (K) = ]im A -+oo K/A. 



These morphological measures are in general useful to characterize the 
structure of various materials, e.g. foams, gels, membranes, granular and 
porous electrode systems. These measures can be calculated for both de- 
terministic and stochastic geometries, and are related through integral ge- 
ometry to Minkowski junctionals [57[. The ensemble average morpholog- 
ical measures now can be related through above integral of mean curva- 
ture, integral of mean square curvature and integral of Gaussian curvature. 
Hence for a large nanostructured surface, Eqs. [9] and [12] are related through: 
(cq) = Hindoo Cd/A. Since the reciprocal of compact and diffuse layers 
local capacitance is additive in nature we will perform averaging over an en- 
semble of random configurations. The ensemble averaged (inverse) capacity 
is now 



Using binomial expansion of relative deviation in capacity, its simplifies to 
following equation 



(1/Cff) + (1/co) 

l/c H + (l/((c G )+6c G )) 



(16) 



l/c H + l/(c G ) + (5c%)/(c G r + 0((5c%)) 



(17) 
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Eq. [T7] for random morphology electrode can be written as the ensemble 
average of inverse of capacitance as: 



where 1/ (cs) = I/ch + 1/ (cg) and (cs) is mean compact and diffuse layer ca- 
pacitance density for average surface morphology The relative mean square 
fluctuation in capacitance ((3 2 = (Sc G ) / (cg) 2 ) is: 

(3 2 = 7 2 ll (HI A? I [l - (H/A)l D - {QP/A) - (K/A)) (l 2 D /2)] * (19) 

where 7 represent the extent of fluctuations in surface morphology 

Using an assumption of ergodic field [56], Eq. [18] in combination with 
Eqs. [12] and [15] simplifies to the mean EDL capacitance in two components. 
The non-fluctuating contribution ((eg)) from this equation is represented as 



(cs) = (l/c H +l/{c G )) 



-1 



ChCq 


1 - (H/A)l D - {{^1 A) - (K/A)) 




ch + c 


\ - (H/A)l D - ((JP/A) - (K/A) 


) vim] 



(20) 

The contribution from fluctuating component (/3 2 /(cg)) in the capacitance 
which arise due to stochastic nature of morphology is 



(5 2 1 YUH/A) 



(cg) c 



1 - (H/A)l D - ({IP I A) - (K/A)) (11/2) 
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(21) 



Eq. [TS] in combination with Eqs. [2D] and [2T] is generalized average capaci- 
tance density (q) for an arbitrary random geometry. This expression relates 
the capacitance to the surface to the averaged mean and Gaussian curva- 
tures of the random morphology electrodes. Hence, the capacitance density 
(cd) of the complex morphology electrode is dependent on an overall surface 
morphological and topological characteristics, viz, H, H 2 and K. 

4-1. Porous Membrane with Fluctuating Pore Sizes 

The complex nature of the nanoporous electrode may results in various 
distributions in surface shapes (curvatures) and pore lengths. In this section, 
we model the influence of intrapore and interpores size fluctuations. Interpore 
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size fluctuations arise due to the presence of polydispersity in pore sizes, 
which is often much larger than intrapore fluctuations. Hence, the nature 
of pore size and shape fluctuations is critical in determining capacitance of 
porous electrodes. For narrowly distributed pore sizes, e.g., carbide derived 
carbons (CDC) micropores with the mean pore size is less than 2 nm, only 
the intrapore local surface fluctuations are significant. There can be low 
intra-pore fluctuations possible without choking pore in the electrode. For 
broadly distributed mesopores not only the fluctuation in interpores sizes 
is important but also the intrapore curvature fluctuations and hence there 
could be simultaneous contributions from both. The porous materials have 
pores of various sizes and shapes, e.g. activated carbons, such systems are 
approximated as membranes electrode with micropores and mesopores. 

Here, we confine our analysis to simple case of nanoporous materials as 
a porous membrane model which has an array of cylindrical pores. These 
pores are allowed to fluctuate weakly, so that the topology of the membrane 
electrode is not changed. Typical geometrical features of the membrane are 
characterized by pore separation distance (w), membrane thickness (pore 
length) (I), the outer macroscopic geometric area of membrane electrode 
(2Aq) and number of pores (Nh). In order to calculate Nh, we assume a 
simple model of membrane with a hexagonal array of cylindrical pores whose 
number of pore is given by Nh = 2A /(y/3(2r + w) 2 )j58|, where r is the 
radius of cylindrical pore. Since the number of pore per unit area is usually 
very large, in orders of 10 12 , the area of two planar surfaces, 2(Aq — irr 2 Nh) 
, connecting two ends of pores of the membrane is small and hence can be 
neglected. The major part of area contribution comes from the area inside 
pores. The total area of such a porous surface is given by area of individual 
pore times pore density and, which is A « 2irrl Nh- For a cylindrical pore, 
mean curvature H = l/2r is constant and the average mean curvature is 
H/A = J HdS/A = l/2r. Similarly, square of mean curvature is: H 2 = 
l/4r 2 and the average square of mean curvature is: H 2 /A = J H 2 dS/A = 
l/4r 2 . The average Gaussian curvature of a porous membrane is calculated 
using the theorem of Gauss-Bonnet as K/A = J KdS/A = 47r(iV c — Nh), 
where N c is number of connected surfaces. For a single membrane electrode, 
the connected component N c = 1 and K/A = (2/rl)(l — Nh)/Nh ~ — 2/r/(for 
large Nh electrode). 

Knowledge of various morphological measures for a given model allows 
us to illustrate the use of Eq. [18] or Eq. [20] under smooth surface limit 7—7-0 
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hence taking ion with the effective ion size in the HL or thickness of HL, r# 
and the diffused layer of the EDL is assumed to start after the Helmholtz 
layer. The effective or attuned pore radii for diffused layer are accounted by 
subtracting HL thickness r# from pore radius r, i.e. r — r#. Substituting 
the morphological quantities in Eq. [201 an explicit equation for the smooth 
pore membrane is obtained as 

, > c H c [1 - l D /2(r - r H ) - l 2 D /8(r - r H f - l 2 D /l(r - r H ) + l 2 D /N h l(r - r H )] 

{CS) c H + c [1 - l D /2(r - t h ) - &/8(r - r H ) 2 - l 2 D /l(r - r H ) + l 2 D /N h l(r - r H )\ 

(22) 

Similarly, the relative mean square fluctuation in capacitance (j3 2 ) due to 
curvature fluctuations in pore sizes or roughness along the contour of pores 
is 

P 2 _ 1 7 2 ^/4(r-r„) 2 



(cg) Co [1 - l D /2(r - r H ) - l 2 D /8(r - r H ) 2 - l 2 D /l{r - r H ) + l 2 D /N h l{r - r H )] ; 

(23) 

It may be noted from Eqs. [22] and [23] for two components of capaci- 
tance density equation, q = ^-^y + ^yj , has a contribution from the 
Helmholtz layer, the diffuse layer and morphological parameters, viz, pore 
radius (r), pore length (/), number of pores (N h ) and coefficient of fluctu- 
ation (7). The capacitance of porous surfaces has a form of rational poly- 



nomial, viz, Pade approximant [59(, in inverse attuned pore size (r — r#). 
Eq. 1221 shows that the ionic concentration dependent Debye length and 
attuned pore size (r — r#) are length scales that affect the capacitance of the 
porous membrane. As mentioned, earlier that the diffuse layer contributions 
in Eqs [22] and [221 have solution property dependent term, geometry depen- 
dent term and a coupling term between solution and geometry. For pore size 
r — th = lnm, the relative importance of these three terms for Id ~ 1 nm 
(for 0.1M electrolyte solution) and is 1:1/2:1/8 and for l D fa 3 nm (for 0.01M 
electrolyte solution) its 1:3/2:9/8. Hence, we conclude that diffuse EDL will 
have a significant influence on the capacitance of a porous system. 

In Figs. 3(a) and 3(b), we analyse the effect of Helmholtz layer thickness 
and the effect of change of concentration (or Gouy length) in the absence 
of the fluctuations in contour of the pore (7 — > 0), respectively. Fig. 3(a) 
shows the attuned (linearized) Gouy-Chapman-Stern (AGCS) capacitance 
density plots for the porous membrane electrode, but without pore diameter 
fluctuations (7 = 0). The curves show the effect of HL-thickness on overall 
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Figure 3: Capacitance plots for a porous membrane electrode with cylindrical pores ob- 
tained from Eq. [22]for 1:1 organic electrolyte with e = 38 at 298 K. (a) Effect of compact 
layer thickness, rjj, on the capacitance profile. The plots are generated using concentra- 
tion, c = 1 M, ejj = 6 and rn = 0.35, 0.4 and 0.5 nm. (b) Effect of concentration on 
capacitance. The plots are generated using rjj = 0.35 nm, ejy = 6 while the concentration 
used in various curves are: 0.01, 0.1 and 1 M. 



capacitance. These plots are obtained using Eq. [22] which show two regions 
in capacitance behaviour and a crossover between them. The point parti- 
cle nature of ions in diffused layer allows us to extrapolate our capacitance 
vs pore-size plot in the narrow pore size regime, but this continuum model 
for the electrical double layer capacitance is valid upto pore size 2r H + £ D 
hence prediction below this value will require dominant contribution from 
the electronic capacitance of the electrode. These observations usually indi- 
cated in several experiments but not captured by heuristic modelsJH, O, 17 



The capacitance becomes independent of the pore size beyond 2 nm, while 
the experiments show dependence up to 10 nm size. This could be due to 
not inclusion of roughness or morphological fluctuations in porous electrodes. 
The quantitative predictions near the highest capacitance region may require 
electronic space charge correction from the electrode as well as other contri- 
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butions like adsorption, excluded volume, etc. in the high field compact layer 
region. 

Fig. 3(b) shows the effect of concentration and pore sizes on the GCS 
capacitance density. Three curves are obtained for different concentration of 
electrolytes (or Id). There is a gradual disappearance and a simultaneous 
shift in the location of the minimum point with the decrease in concen- 
tration. Other conclusions can be drawn from these curves are: (a) high 
concentration result, Kr > 1, has a minimum value of capacitance at pore 
size 2r = 2th + ft -1 , (b) low concentration results, kt < 1, there is gradual 
reduction of capacitance. Anomalous increase in capacitance is observed in 
CDC where pore size is less than 1 nm [14j but this cannot be explained 
purely on the basis of ionic capacitance. The capacitance has a sudden rising 
region followed by a gently rising with increasing pore size and a crossover 
between these two regions is seen here. This is an interesting finding for 
system with unimodal pore size and influence of the fluctuation in pore sizes 
in capacitance behavior will be shown in Fig. 4. 

In order to see the effect of morphological fluctuation on capacitance den- 
sity, we plot capacitance of a porous membrane electrode. The pore shape 
and size can fluctuate along the contour of the pore. Fig. 4 shows the effect of 
surface morphological fluctuations on the capacitance of a porous membrane 
electrode. Fig. 4(a) describes capacitance vs mean pore size in the presence 
of small fluctuations in the micropore regime. We observe from the curves 
that as we increase the value of surface morphological fluctuation parameter 
7, the sudden rise region becomes weaker while the slow rise region becomes 
slightly faster rise. All the curves merge around 2 nm and no effect of small 
fluctuation is seen beyond this. Thus, the small fluctuations in micropores 
are important upto 2 nm. Fig. 4(b) shows the effect on capacitance due to 
large fluctuations in mesopores regime. As we increase 7, the value of capac- 
itance is lowered, compare to smooth porous electrode. The effect of large 
fluctuation is seen for pore diameter much larger than 2r#+ Id- Increasing 
the mean pore size in mesopore regime reduces the influence of fluctuation, 
and finally, all the curves merge around 10 nm. No curvature effect is seen 
beyond 10 nm, which suggest that large morphological fluctuation enhance 
capacitance dependence on mean pore size upto 10 nm. This is in agree- 
ment with several experimental reports of capacitance vs pore size data (see 
comparison with experimental data). The small fluctuation (7 < 2) affects 
the micropore regime while for large fluctuations (7 > 2), affects the capac- 
itance in the mesopore regime. For pore diameter less than 2 nm, the slope 
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Figure 4: Effect of the morphological fluctuations in a porous membrane system on ca- 
pacitance vs mean pore size plots, (a) Small fluctuation in pores sizes with 7 = 0, 0.3, 
0.6 and 0.9. (b) Large fluctuation in pores sizes with 7 = 2, 4, 6 and 8. All plots are 
generated using r# — 0.35 nm, en = 6, for 1:1 organic electrolyte with e = 38 at 298 K. 
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of the rapid rise region is almost same. In case of mesopores the fluctuation 
results in a decrease of capacitance. Larger is the fluctuation, larger is the 
pore size up to which the effect is seen. The capacitance increases rapidly 
in the transition region between micro and mesopores followed by a gradual 
increase in mesoporous regime and limiting plateau value of a flat electrode 
is observed at large mesopores. Hence large morphological fluctuations play 
an important role to enhance the capacitance density dependence on mean 
pore size. 

4-2. Forest of Nanorods with Morphological Fluctuations 

In this section, we obtain capacitance for the forest of nanorods and com- 
pare them to that of nanopores. These two systems differ only in convexity 
of their electrode surface. The typical dimension of the forest of nanorods 
electrode is characterized by rod separation distance, (it;), electrode thickness 
(rod length), (/), the geometric area (A ) and number of rods (N r ). In order 
to calculate N r we assume an electrode with a hexagonal array of cylindrical 
rods whose number of rods is given by JV r = 2A /(v / 3(2r + w) 2 ), where r 
is the radius of cylindrical rod. The total area of such a forest of nanorods 
structured surface is given by area of individual rod times rod density (N r ) 
and it is A « 2nrl N r . Our theory is applied to nanoforest electrode with in- 
ter rod separation (w) which can accommodate compact as well as a diffused 
layer between them, viz, w 3> 2r# + Id- The capacitance for array of cylin- 
drical rods is obtained by replacing, H/A = (— l/2r), H 2 /A = (— l/2r) 2 , 
and attuned rod size for diffused layer is obtained by adding compact layer 
thickness to rod radius: r + th- Here the diffuse layer starts after a dis- 
tance r + th from the center of each rod. The average Gaussian curvature 
of the nanorods forest is calculated using the theorem of Gauss-Bonnet as 
K/A = j KdS/A = Atc(N c - N h ) /A. Using Eq. [20]we obtained the mean ca- 
pacitance for N r electrically connected forest of nanorods: N c = N r , Nh = 
and K/A = (2/rl). 

, v chop [1 + l D /2(r + th) - l 2 D /8(r + r H f + l 2 D /l(r + r H )} 

[CS)R c H + Co[l + l D /2(r + r H )~l 2 D /8(r + r H ) 2 + l 2 D /l(r + r H )] { ] 

The relative mean square fluctuation in capacitance which arise due to surface 
fluctuations over rods is obtained using Eq. [21] as 

/? 2 _ 1 7 2 Wr + r H ) 2 



(c G ) c [1 + l D /2(r + r H ) - l 2 D /8(r + r H ) 2 + l 2 D /l(r + r H )] 
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Fig. 5(a) shows the effect of HL thickness on the GCS capacitance density 
vs rod diameter plots in a forest of nanorods electrode where 7 = 0. Plots are 
made using Eq. [2H As we increase the HL thickness, the value of capacitance 
decreases. The HL increases the attuned radii curvature of rods and hence 
resultant decrease in capacitance contributions from the diffuse layer. Fig. 
5(b) shows the effect of diffuse layer thickness, which increases with decrease 
in the concentration of electrolyte, leads to decrease in capacitance. Thus 
for nano rods with the larger attuned radii of curvature, effect will largely be 
seen in the capacitance for dilute solution. In case of concentrated solution, 
it simply increases the capacitance but no enhancement in capacitance is 
seen in microrods. 

Fig. 5 (b) shows the effect of concentration on the capacitance of as- 
sembly of nanorods electrode using Eq. [2H (7 = case). Unlike cylindrical 
pore geometry, the capacitance of rod geometry is enhanced with reduction 
of diameter. The overall capacitance is very much like the Helmholtz ca- 
pacitance, but it is enhanced mainly for small rod diameter. Increase in 
concentration results in higher capacitance. We compare the capacitance ob- 
tained in porous membrane electrode with nanorods with same surface area, 
but different (local) convexity can have qualitatively different capacitance 
dependence on size. These are purely geometric effects of the surface curva- 
tures. Thus, it can be concluded that membrane electrode with cylindrical 
nanopores have different capacitance as compared to forest of cylindrical 
nanorods although both have same geometrical surface area. We see from 
the plots that the curvature strongly affects the diffuse layer capacitance. 
Similar studies by Compton and coworkers 32|, |52|] has shown that for hemi- 
spherical and nanotubes electrodes, there is enhanced diffuse double layer 
capacitance due to increase in curvature. In fact, our analytical model cal- 
culations also show such an increase in diffuse double layer capacitance. But 
the introduction of compact (Helmholtz) layer capacitance in series with the 
diffuse layer capacitance reduces the total capacitance. Here overall capac- 
itance is controlled by smaller of the two. It is important to note that the 
curvature effect is very much reduced in case of nanorods electrode unlike 
the diffuse layer capacitance without attuned radii correction. Furthermore, 
as we decrease the concentration and hence increase in the Gouy length, we 
see the curvature effect to capacitance reduced. The diffuse layer is extended 
in low concentration regime and may require accounting for overlapping of 
EDL @,[6l|. 



Fig. 6 shows the effect of morphological fluctuations on the capacitance 
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of forest of nanorods. Increasing the curvature fluctuation parameter (7). 
The capacitance behaviour of nanorods changes qualitatively from smooth 
nanorods capacitance behavior to nanoporous capacitance behaviour. For 7 
= 2 , the curve nearly approaches a horizontal line of a planar electrode, and 
effectively no curvature effects are seen. Curves show lowering in capacitance 
with an increase in the value 7 and curves appear similar to porous electrode 
capacitance with a gradual fall in smaller rod thickness. The physical reason 
is, for an assembly of nanorods electrode with large curvature fluctuations on 
the nanorods results in formation of pits or pore like structure which made 
the capacitive response to behave like a porous system. 



5. Comparison of Theoretical Results with Experimental Data 

Fig. 7 shows the theoretical and experimental capacitance behavior of 
organic electrolyte TEABF4 in CH3CN in different pore size regimes and 
concentrations. Data were taken from references [H, Q, 17]. Dielectric con- 



stant of CH3CN is taken to be 38 at 298K, bare ion diameter of TEA + is 
0.68 nm, BF 4 ~ is 0.46 nm and solvent diameter is 0.34 niii|02 . We use in 
our calculations for the compact layer dielectric constant e# = 5.3 and the 
bare ion radius of TEA + as the compact layer thickness r# = 0.34 nm. The 
plots are obtained using Eq. [20] for smooth pores and Eq. [18] in combination 
with Eqs. [22] and [23] We have mentioned earlier that we use bare ion size 
for the compact layer thickness and a reasonable value for the compact layer 
dielectric constant to evaluate four theoretical curves, as an illustration, of 
overall behavior of capacitance predicted by the theory and compared them 
with capacitance behavior in some of the recent data. The theoretical plots 
for capacitance clearly indicate the minimum in an organic electrolyte. Black 
line is for 1.5 M electrolyte and electrode without morphological disorder, i.e. 
7 = and black broken line are for same concentration electrolyte but the 
electrode has some morphological disorder, i.e. 7 = 2. Similarly, two lower 
curves have the electrode with same morphological disorder parameter 7 = 6 
but differ in electrolyte concentration. Two electrolyte concentrations 1.4 M 
and 1M, are compared with experimental data to illustrate the point that 
capacitance is dependent on electrolyte concentration. TEABF 4 in CH 3 CN 
at 1.4 M have higher capacitance than 1 M TEABF 4 in CH 3 CN for porous 
electrode with 7 = 6. In agreement with the theoretical prediction, the effect 
of morphological fluctuation is seen upto 10 nm in the experimental data. 
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Experimental comparison can further improve with inclusion of the con 



tribution from electronic space charge capacitance of the solid 63|, [64], [65 



as 



its contribution will be important in understanding the anomalous rise region 
and its maximum, but we have ignored this in the present calculations. 



6. Summary and Conclusions 

In summary, a model for EDL capacitance in curved nanostructured 
electrode is proposed within the frame work of classical (linearized) Gouy- 
Chapman-Stern theory. Our analytical theory for electric double layer ca- 
pacitance of nanostructured and porous electrodes of arbitrary geometry is 
applicable to the wider class of electrode morphology and electrolyte sys- 
tems. Theory for arbitrary electrode geometry is developed using the "mul- 
tiple scattering formalism" in surface curvature, which includes geometrical 
as well as topological information about the surface. The region near the 
electrode is divided into a strong electric field region- the compact layer and 
the region beyond the compact layer with a relatively weak electric field 
region- the diffuse layer (assumed to obey linearized Gouy-Chapman equa- 
tion). The capacitance in the diffuse layer is expressed as a function through 
convergent expansion in power of the ratio of Debye length to the principal 
radii of curvatures of the interface. The contribution of compact layer and its 
size correction in attuned geometric size in diffuse layer has a strong effect on 
the local capacitance. This is introduced through the capacitance in series 
behavior of local compact layer and local diffuse layer contributions. Theory 
is further extended for realistic electrodes with a ubiquitous morphological 
disorder /roughness which is characterized through variance in mean curva- 
tures. Our theory is general enough to handle arbitrary morphology of an 
electrode, but we illustrate its significance for the simple system of forest of 
cylindrical nanopillars and nanoporous membrane with added complexity of 
morphological disorder. Theory allows us to analyze the effect of concen- 
tration, compact layer thickness, curvatures and its fluctuations along the 
pore. 

This paper further establishes for the capacitance of nanostructured elec- 
trodes the following points. The theory predicts a minimum value of the 
capacitance at minimum allowed mean pore size 2rn + Id an d a gentler rise 
towards larger pore sizes with a limiting capacitance plateau (at planar) as 
observed in experiments. The capacitances for nanoporous membrane elec- 
trode and nano-forest electrode, with same geometric area, have qualitatively 
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different capacitance dependence on sizes. Hence, the overall convexity of 
electrode surface has a profound effect on the ionic capacitance. The cur- 
vature fluctuations in nanostructured electrodes have profound influence on 
the capacitance in the mesoscopic region. The small curvature fluctuations 
along the contour of pore or intrapore geometric fluctuations affect the rise 
of capacitance which extends upto 2 nm size. Hence, weak intrapore fluctu- 
ations slightly enhance the pore size dependence of capacitance. The effect 
of large curvature fluctuations manifest themselves through prolonged de- 
pendence of capacitance on pore sizes and it persist in up-to nearly 10 nm 
in mesopore regime. Surprisingly, the large fluctuations in nano-rods forest 
electrode results in qualitatively different dependence on rod thickness. 

Finally, our model is a step forward to develop a general theory of capaci- 
tance in complex interfacial systems. This theory can be further extended by 
accounting quantum-mechanical contributions from electronic space charge 
in the electrode as well as more details in description of ion-ion correlations 
of EDL. These extensions of the theory may help to understand, the yet un- 
solved anomalous increase in capacitance in pore size less than 1 nm so called 
the "Gogotsi-Simon effect". Analogous problems of finding effect of non- 
ideality of electrode surface, compact layer disorder and dynamic frequency- 
dependent response of nanostructured and porous electrodes will be reported 
elsewhere. 
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Appendix 

The Debye-Hiickel or linearized Gouy-Chapman equation for the potential 
relative to the bulk solution <p(r) is: 

(V 2 - K 2 )<f){r) = 0. (A. 1) 

For a domain D bounded by a conducting surface S which is held at 
constant potential <p + \s = 4>o and far away from surface , viz. <p(r — > oo) = 
0. The Green's function (GF) G for a linearized PB Eq. IA. II satisfies 

(V 2 - K 2 )G(r, r') = -5 3 (r - r') (A. 2) 

with the homogeneous boundary conditions at surface S and far away from 
surface (r — > oo), viz. G\$ = 0. The solution of Eq. IA. 21 for the potential 



relative to the bulk solution, <f)(r) can be formally written as 43 

dG(r,p+) 



dSn 



drif 



o / d 3 r'V 2 G(r,r') 
Jv 

-k 2 [ d 3 r'G(r,r') ) (A. 3) 



v 



We seek to express G in terms of Go(r,r') = (l/47r|r — r / |)e _K ' r_r ''', the 
free space Green's function satisfying Eq. IA. 2l for the entire infinite domain. 
The Green's function G is expressed as |43| . 



G(r,r') = G (r,r')-2 J G (r, 



dG (a,r') 
— dn 



+2 J G {r,a) — — — dS a dS p 

+ ■•• (A. 4) 



Using Gauss's law and the fundamental singularity at the boundary [42|, |43 
the differential capacitance density is 

da) = — / d A r 



47r ./,/ dn r 
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^ / dV 



dGn(a, r') 



-2 



dn a 

dG (a,/3)dG (^r') 



dS. 



P 



+2 



dn a dnp 
, , dG (a,P)dG ((3,~f)dG (~f,r>) 



dn a 



dn 



P 



dn 



dSpdS^ 



7 



(A. 5) 



The terms in Eq. IA. 51 can be looked upon as one-, two- and three-scattering 
terms as 



Si 



, 3 ,dG (a,r') 



d A r 



dn a 

dG (a,P)dG (P,r') 
dn a dnp 

dGo(a,/3)<9G (/3,7) dGo^r') 



dS B d 3 r' 



dn a 



dnr- 



dn~. 



dSpdS y d 3 r' 



(A. 6) 
(A. 7) 
(A. 8) 



Eq. IA. 61 may be expanded through the local surface coordinates a, (3, 7. For 
a weakly curved surface where Debye-Hiickel screening length is much smaller 
than the smallest scale of curvature the scattering kernel dG ((3 + ,a)/dn a is 
expressed through a local coordinate system 42j, |43( with the z-axis parallel 
to inward normal vector n a and a tangent plane on which projection is made. 
The local equation of surface, S in terms of curvature radii R\(a) and ^(a) 

z a = (1/2) (x 2 /R 1 (a) + y 2 /R 2 (a)) +■■■ (A. 9) 

is introduced to Eq. IA. 6l through surface area element dSp = y/gdxdy where 
g = 1 + (Vz a (x, y)) 2 . Using d/dn a = —d/dz the kernel dG (a, (3)/dn a under 
planar approximation [42| reads to first order as (dG (a, (3)/dn a ) = —{zj 
p)(dG (p)/dp), where p = \a' — f3'\ = (x 2 + y 2 ) 1 ^ 2 is distance in tangent plane 
, Go = exp(—p)/4Trp the Green's function in tangent plane. Now we can 
rewrite the one-scattering integral Eq. IA. 61 as 



Si 



dxdy 



r lz 2 dG 



2 p dp 

which is further simplified using the angular averages (< 



> 



\Z a ) 



2n 



dB 



cos 2 9 sin 2 9 
+ 



7T 

2 



1 



RAa) 



Ri(a) 



(A. 10) 



(A. 11) 
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1 , , x 1 f 2n .„ ( cos 2 6 sin 2 6 



P^ Zl) ~ Jo d9 \R 1 (a) + R 2 (a) 

3?T / 1 1 \ TT 

T V^i(a) 2 + ^2(a) 2 J + 4i?i(a)i2 2 (a) ( ' ^ 

Substituting Eq. IA. lJ to lA. lOl and integrating over p we finally get the one- 
scattering term as 

Sl = k ~ ^ ( 2j R( a ) 2 ~ 2 J Ri(a) J R 2 (a)) + ° (#V ^' 13) 

where F Q = l/fl(a) = (l/2)(l/ J R 1 (a) + l/R*(a)) and # Q = l/ J Ri(a) J R 2 (tt). 
Similarly on iteration in Eq. IA. 71 and Eq. IA. 81 two- and third-scattering 
terms are obtained as 

o(-L), Es=ii^T» + of4 > ) ( A - 14 ) 



« 2 fl(a) \R 3 J ' k 3 R(a) 2 \R 3 / 

Now substituting Eq. IA. 131 and Eq. IA. 141 to Eq. IA. 5| after simplification 
the capacitance density at position a is 

1 (A. 15) 



c(H a , K a ) = — 

47T 



1- 1 -rf 
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Figure 5: Capacitance plots for an assembly of cylindrical rods with no morphological 
fluctuation (7 = 0) obtained from Eq. [22] Fig. (a) shows the influence of the compact 
layer thickness (ru), e.g. 0.34, 0.4 and 0.7 nm. Fig. (b) shows the influence of the 
electrolyte molar concentration, e.g. 0.01, 0.1 and 1 M. All the plots are generated using 
th = 0.34 nm. en = 6 for 1:1 organic electrolyte with e = 38 at 298 K. 
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Figure 6: Effect of morphological fluctuations on the capacitance of a nanorods forest 
electrode. Here fluctuation parameter is varied as: 7 = 0.3, 0.6, 0.9, 2, 4, 6, 8. All the 
plots are generated using r# = 0.25 nm, e = 6, for 0.05 M H2SO4 with e = 78.6 at 298 K. 
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Figure 7: Theoretical capacitance curves for a porous membrane electrode model of the 
electrode with and without morphological disorder. Experimental data points are from ref 
[3 51 Hi- ( a ) Black line is for 1.5 M tetraethylammonium tetrafiuoroborate (TEABF 4 ) 
organic electrolyte in acetonitrile (ACN) with e = 38 at 298 K, e# — 5.3, rn = 0.34nm, 
7 =0. (b) Broken black line is the same case as in (a) but with different 7 — 2. (c) Red 
line is for 1.4 M concentration and 7 = 6 while the rests of parameters are same as in (a), 
(d) Blue line for 1 M concentration and rest of system parameters are same as in (c). 
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